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ABSTRACT 

We compute numerically the scalar- and tensor-mode induced Stokes 
parameters of the cosmic microwave background, by taking into account the 
basis rotation effects. It is found that the tensor contribution to the polarization 
power spectrum get enhanced and dominates over the scalar contribution for 
low multipoles in a universe with or without recombination. Furthermore, 
we show that all full-sky averaged two-point cross-correlation functions of the 
Stokes parameters vanish. We thus comment on the cross-correlation between 
the anisotropy and polarization of the CMB, and calculate the expected signal 
to noise ratio for the polarization experiment underway. 

Subject headings: cosmic microwave background — cosmology: theory — 
polarization 
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1. Introduction 

The detection of the large-angle anisotropy of the cosmic microwave background by the 
COBE DMR (^moot et al. 199^ ) is an important evidence for the large-scale space-time 



inhomogenities. Since then, several small-scale anisotropy measurements have hinted that 
the Doppler peak resulted from acoustic oscillations of the baryon-photon plasma on the 
last scattering surface seems to be present. It is believed that measuring CMB anisotropics 
on all angular scales can give invaluable information of the early universe, such as testing 
inflation, pinning down cosmological parameters, and differentiating scalar from tensor 
perturbations (see, e.g., [Steinhardt 1994|) . 



Complementary to the anisotropy, polarization is another important property of the 
CMB. The degree of CMB polarization has been extensively calculated (see |[\lg fc Mg 1996 



and references therein). In a universe with standard recombination, the polarization is 
about 10% of the anisotropy on small angular scales while the large-scale polarization is 
insignificant. If early matter reionization of the universe occurs, the large-scale polarization 
would be greatly enhanced to a 10% level with the large-scale anisotropy reasonably 
unaffected. Meanwhile, the small-scale anisotropy would be significantly suppressed. 
Thus, measuring CMB polarization might provide additional information about the early 
universe. So far, only experimental upper limits have been set ( |Lubin, Melese, fc Smoot" 



[Partridge, Mowakowski, fc Martin 1988| ; [VVollack et al. 1993| ), and the signal level is 



expected to be achieved through the use of new technology and instrument design (P. T. 
Timbie 1995, private communication). 

In previous calculations of the degree of CMB polarization, approximations have more 
or less been made. For instances, the separation angle in the two-point auto-correlation 
functions has been assumed small, or/and the basis rotation from the k-basis of the 
perturbation to the fixed basis of an experiment has been totally neglected. Although it is 
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expected that these approximations would be vahd for small-angle polarization calculations, 
they may break down on large-scale calculations. In this Letter, we will compute the Stokes 
parameters by taking into account the basis rotation effects. Also, we will consider the 
two-point cross-correlation functions of the Stokes parameters. 



Stokes Parameters of the CMB 



Polarized light is conventionally described in terms of four Stokes parameters: 
{I,Q,U,V), each of which is a function of the observational position in the universe, r, 
as well as a unit vector pointing to the celestial sphere, e. In polarization measurements, 
it is convenient to choose the North Celestial Pole as the reference axis, z, and the linear 
polarization of the CMB is measured in terms of Q and U defined by: Q = T/v,5 — Te^w and 
U = Tne^sw ~ Tnw,se, where T/v,5 is the antenna temperature of radiation polarized along 
the north-south direction and so on ( [Lubin, Melese, fc Smoot 1983 ). Let us define T = I — I 
as the temperature fluctuation about the mean. It is known that V would decouple from 
the other components. Thus, we will calculate the Stokes components (T, Q,U) of the CMB 
induced by scalar and tensor perturbations. We first decompose (T, Q, U) into Fourier 
components (T^, Qkj f^k)? which are induced by a Fourier mode with wavevector k. Then, 
for a fixed k, we solve the radiative transfer equation for the Stokes parameters (T^, Q^; ^k); 
which are defined with respect to the k axis. The primed and unprimed parameters are 
related by a simple rotation ( Bond fc Efstathiou 1987 ), 




cos 2il) sin 2il) 
— sin 2il) cos 2%l) 

and Tk = T^, where the rotation angle if) is given by 

k • z — (k • e)(z • e) 
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cos iIj 



1 - fk ■ e) 
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1 - fk ■ e) 



(2) 



z- e 



- 5 - 



In the tensor-mode case, = (l/2)a(l — /i^) cos 20, Q'^. = (l/2)/5(l + ;U^)cos20, and 
f/^ = /?/isin20 for the + mode solution, where (/i = cos^,0) are the polar angles of the 
unit vector e in a coordinate system with k as the polar axis, a and (3 are functions of fi 
only, usually expanded in terms of Legendre polynomials: a{fi) = I]/(2^ + l)aiPi{fi) and 
/3(yu) = X]i(2^ + ^)PiPi{fi'), where a/ and Pi have previously been calculated (see [Ng fc Ng 
|1996| and references therein). The x mode solution is given by the same expressions except 
replacing cos 20 by sin 20, and sin 20 by — cos 20. In the scalar-mode case, = a, Q'y_ = (3, 
and t/^ = 0. 



3. Anisotropy and Polarization Power Spectra 

In CMB measurements, the measured Stokes parameters are usually expanded in terms 
of the spherical harmonics, 

T(r,e) = J2^irnir)Yimie), Q{r,e) = J2 1irnir)^imie) , f/(r,e) = ^ Mi^(r)Y;™(e). (3) 

l,m l,m l,m 

From this, we define the anisotropy power spectrum, Ci = Ti, and the polarization power 
spectrum, Cf = Qi + Ui, where 

2/ + 1 r 

Ti = {J2tLtim) = J de^de^Piie, ■ 62)^X2), (4) 

and Qi and Ui are given by similar expressions. In equation (^, the angular brackets denote 
taking an average over all observational positions, i.e., (/(r)) = l/(27r)'^/ drf{r), and Ti 
represents T(r,ei). 

In the following, we will present the results for the anisotropy power spectrum, 
Ci = Ci-^"^ + Cf^\ and the polarization power spectrum, C[ = Cf^^^ + Cf^^\ where S 
and T represent respectively the scalar and tensor contributions. Specifically, we work in a 
flat CDM model with no recombination, in which the Hubble parameter h = 0.5 and the 
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baryon density fi^/i^ = 0.0125. Also, we assume that Ca^^) = C2^'^\ the scalar power index 
is Hs = 0.85, and the tensor's is rit = —0.15. 



To evaluate the power spectra in equation (g), we have used the Boltzmann code 
fc Ng 1996| ) to solve for ai and For the anisotropy, the integral Ti can be easily evaluated 



and analytic solutions for Ci^^^ and C/"^^ can be found in terms of a^'s ( [Crittenden et al] 
1993|) . For the polarization, unfortunately, the multi-dimensional integrations involved in Qi 



and Ui are difficult. However, it has been attempted to calculate Qi and Ui approximately. 
In both scalar and tensor cases, by neglecting the basis rotation, i.e., setting ip = 0, 



analytic expressions for Cf and Cf have been obtained by |Mg fc Mg 1995| . Using these 



approximate formulae, we have plotted in Figure 1 the scalar contribution, Cf /Ci, the 
tensor contribution, Cf /Ci, and the total polarization. 

To study the effect of the basis rotation, we numerically compute the integrals Qi 
and Ui using the Gaussian quadratures integration method. We have computed the power 
spectra for / = 2, 3, 5, 10, 20, 30 with a 10% accuracy (improving the accuracy by an 
order of magnitude would require much more computer time). We have found that the 
coordinate-dependent Stokes parameters Qi and Ui are equal for a fixed /, as expected for 
Gaussian random perturbations. For the scalar mode, including the rotation angle does 
not apparently change the results from using the approximate formula with ip = 0. On 
the contrary, in the tensor case, the multipoles for / < 30 are modified, and the above 
approximation is valid only for / > 30. This is shown in Figure 1 with solid symbols. 
Note that Cf is enhanced by an order of magnitude and now exceeds Cf , and the 
tensor contribution is dominant for I < 4. This is contrary to the result in Figure 3 of 
Crittenden, Davis, & Steinhardt 19931 where the tensor contribution is subdominant for all 



I. On the other hand, for a universe with standard recombination, the tensor contribution 
is predominant for multiples with / < 30 ([Crittenden, Davis, fc Steinhardt 1993 ; Ng fc 
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Ng 1996). Thus, the tensor mode plays an important role in the large-scale polarization. 
Although the basis rotation effect raises considerably the low multipoles of the tensor 
spectrum, it only slightly alters the r.m.s. degrees of polarization that have been calculated 
in |Ng fc Ng 199"^ . 

Now, we can construct the auto-correlation functions. Spatial isotropy guarantees that 

iKrJl'm') = 2/ _^ ^^ll'^rnrn'- (5) 

Hence, the temperature auto-correlation function is given by C(0) = (T1T2) = 
(l/47r) J2iCiPi{cosQ), where G is the separation angle. For the polarization, both Q and U 
only have cylindrical symmetry, so qim and Uim do not have a relation like equation (|^). As 
such, we first take a whole-sky average which is denoted by the curly brackets, and then an 
average over all observational positions, namely, 

c^(e) = ({Q1Q2 + = -^^cfPKcose), (6) 

where we have used {y5^(ei)Yj/m'(e2)} = l/(47r)Pi(cos 0)5n"^mm'- In actual observations, 
Yim in equation (|^) should be replaced by WimYirm where Wim is a window function 
appropriate to a particular experiment. However, for many polarization measurements the 
window function only depends on / and is simply inserted in equation (|]). 



4. Cross-correlation Functions 

It is interesting and useful to know the cross correlations among T, Q, and U both 
theoretically and experimentally. The two-point cross-correlation functions (T(r, e)Q(r, z)), 
(T(r, e)?7(r, z)), and ((^(r, e)f/(r, z)) have been considered, where the axes specifically 
chosen to define the Stokes parameters in the z direction are x and y ( Coulson" 



Crittenden, fc Turok 1994 prittenden, Coulson, fc Turok 1995| ; [Kosowsky 1996|) . These 
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cross-correlation functions with this particular choice of axes have limited use in the sense 
of actual observations. Essentially, this is because a single measurement of the Stokes 
parameters at two specific directions e and z has a large sample variance. However, this 
variance can be reduced by surveying a larger part of the sky, and it approaches the 
inevitable cosmic variance for a whole-sky coverage. As such, it is valuable to consider the 
whole-sky averaged {{T1Q2}), ({rif/2}), and {{QiU2}). 

Let us consider {{T1U2}) first, from equations (||) and it is given by 

{{T1U2]) = ^T.iT.tLnim)Pi{cose), 

(EC«/r.) = ^/rfMe2P/(ei-e2)m[/2). (7) 

For the scalar mode, we obtain 

21 4- \ r r 

(EC««m) = ^E(2^i + 1)(2^2 + 1) / k^dk ai.Pi, / de^de2dk 

Piiei ■ e^)PUei ■ k)P/,(e2 ■ k) sin2^2. (8) 

Using the properties that sm2ilj is odd under the parity transformation: e — e, and even 
under the transformation: k — k, one can easily show that {J^m.tim'^im) in equation (||) 
is identically zero, hence, {{T1U2}) = 0. Similar argument shows that {{Q1U2}) = 0, 
but it does not guarantee that {{T1Q2}) vanishes. However, by direct evaluating the 
multi-dimensional integral, we have found that {{T1Q2}) = 0. For the tensor-mode 
case, similarly, using the parity properties of the angles 6, 0, and ip, one can show that 
{{T1U2}) = 0, and {{Q1U2}) = 0. Notice that the interference terms between the cross 
and plus modes in the integrals vanish as expected. For {{T1Q2}), we have evaluated the 
integral numerically and the result is consistent with zero. 

As a result, all full-sky averaged two-point cross-correlation functions vanish. This 
does not mean that there is no correlation between anisotropy and polarization. In fact. 



- 9 - 



the correlated part has been extracted, and its polarization pattern has been correlated 
with the cold and hot spots of the microwave sky ([Coulson, Crittenden, fc Turok 1994| ; 



[Crittenden, Coulson, fc Turok 1995| ). The zero result here only suggests that the correlated 



part has inherent symmetry over the whole sky. It is interesting to estimate how much the 
correlation is suppressed by averaging over an area of sky big enough to reduce the sample 
variance to acceptable levels. However, it can hardly be done without knowing the inherent 
symmetry, and the result should depend on sky sampling. Therefore, it is important to 
find out the exact (no approximations or whole-sky averaging) correlation functions of the 
Stokes parameters between any two arbitary points on the celestial sphere. This would 
reveal the symmetry, and be useful to measurements. 

We have tried to calculate a special case that is (T(r, z)(5(r, e)). For the scalar case, 
the result is 



(r(r,z)g(r,e)) = 2nJ2i2l + 1)0^"^ P^{cosi 



l>2 



(I - 2)! 



dxPi\x)Pv{x) 

1 



k^dka^pi', (9) 



which is similar to the result given in [Coulson, Crittenden, fc Turok 1994| except that 



our result has no azimuthal dependence. This is expected because Q has a cylindrical 
symmetry. Equation (j^) would be useful to the on-going polarization measurement that 
would observe an annulus of 10 spots at constant declination (P. T. Timbie 1995, private 
communication). If the spots are well separated, then the measured anisotropy-polarization 
correlation will be 

ToQ = {TQie))±^, (10) 
/n 



where Q = J27=i Qi/^y and aq are the anisotropy and polarization variances respectively 
To is the measured polar temperature anisotropy and n = 10. (We have neglected the 
detector noise.) In a reionized universe, {TQ{6 = 50°)) ~ — O.OlcrycrQ ( [Coulson, Crittenden. 



Turok 1994[ ), hence the expected signal to noise ratio S/A^ ~ 0.03. When 6 = 4°, 
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S/N ~ 0.3. Whereas, the S/N ratio of the polarization variance (Q^) is -\/5. Note that one 
is hkely to find a relatively large Q when Tq is not at a hot or cold spot on the microwave 
sky. 



5. Conclusions 

We have calculated all whole-sky averaged two-point correlation functions of the 
Stokes parameters of the CMB. Incorporating the basis rotation effects, we have found 
that the tensor contribution to the polarization power spectrum get enhanced for low 
multipoles whereas the scalar contribution remains unchanged, and that all full-sky averaged 
cross-correlation functions vanish. A special case of (TQ) has been calculated, and used to 
show that the signal to noise of (TQ) is low for the polarization measurement underway. 
We intend to calculate the exact correlation functions that would be important to future 
polarization experiments. 
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Fig. 1. — Ratio of polarization multipole to anisotropy multipole vs. I in a standard CDM 
model with no recombination. The short-dashed, long-dashed, and sohd curves denote 
respectively the tensor, scalar, and total contributions, using the approximate formulae with 
no rotation angles. The solid triangles and squares denote respectively the corrected tensor 
and total contributions. 



